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We study causal hydrodynamics (Israel-Stewart theory) of gauge theory plasmas from the 
AdS/CFT duality. Causal hydrodynamics requires new transport coefficients (relaxation times) 
and we compute them for a number of supersymmetric gauge theories including the AT = 4 SYM. 
However, the relaxation times obtained from the "shear mode" do not agree with the ones from the 
"sound mode," which implies that the Israel-Stewart theory is not a sufficient framework to describe 
the gauge theory plasmas. 
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O ■ I. INTRODUCTION AND SUMMARY 

<d ; 

The AdS/CFT duality is a powerful tool to study hydrodynamics of gauge theory plasmas, and it 
has interesting implications even to quark-gluon plasma (QGP). (See Refs. [l|, H, H[ for reviews.) One 
\& ■ robust prediction is a universally small ratio of the shear viscosity r\ to the entropy density s at large 

J> ' 't Hooft coupling: 

\D ' 

g>: s = w (L1) 

. Similarly, one can compute the other transport coefficients such as bulk viscosity, speed of sound, and 

thermal conductivity. 

C . However, standard hydrodynamics (first order formalism) has severe problems such as acausal- 

ity. The first order formalism has the other problems: Equilibrium states are unstable under small 
perturbations Q and the diffusion equation is inconsistent with sum rules [f| . 

One can restore causality, but one is forced to introduce a new set of transport coefficients. Such 
r> | a theory is known as "causal hydrodynamics" or "second order formalism." At present, there is no 

unique formalism for causal hydrodynamics. But probably the most used formalism is the "Israel- 
Stewart theory" E l El- (See Refs. d,Q for reviews.) Another well-known candidate is the "divergence 
type theories" [IOjLUI, which has more attractive features mathematically. In this paper, we focus on 
the Israel-Stewart theory. 

New coefficients which appear in the Israel-Stewart theory may become important in the early stage 
of QGP formation, and in fact it has been widely discussed in the context of heavy-ion collisions. 
For example, a number of groups recently re por ted the results of the (2 + l)-dimcnsional numerical 
simulations of causal hydrodynamics [HI, fl3l. Il4l. HH . Unfortunately, little is known about these 
coefficients: They have been evaluated only for the Boltzmann gas (dilute gas approximation). 

The aim of this paper is to determine these coefficients from the AdS/CFT duality. The AdS/CFT 
duality cannot directly compute these coefficients for QCD, so we compute them for various super- 
symmetric gauge theories including the Af = 4 SYM to see if there is any universality or generic 
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features. Wc determine these coefficients by solving perturbation equations in the Schwarzschild-AdS 
black holes (SAdS) in various dimensions. 
Our results are summarized as follows: 

1. One of the new transport coefficients is tv, which is the relaxation time for the shear viscous 
stress. This coefficient appears both in the "shear mode" and in the "sound mode," but their 
values do not coincide. This suggests that the Israel-Stewart theory is not sufficient to describe 
the gauge theory plasmas. 

2. If one trusts the value of obtained from the sound mode, 1 r n ~ 0.2 fm (for T _1 = I fm.) 
Using the AdS/CFT value of rj/s, this value is not far from the Boltzmann gas estimation. 

3. The numerical values of tv are similar among the theories we consider. 

Explicit results can be found in Sec. II VI In addition, we obtain the relaxation time tj for the charge 
diffusion in those theories. The coefficient r n has been reported in Ref. [ij} for the Af = 4 SYM using 
an expan ding plasma. Wc compare our results and remark implications in an appropriate place (See 
also Ref. [U) 

In the next section, we illustrate the idea of causal hydrodynamics using a simple example. For the 
technical details used in this paper, see App. |XJ We set up perturbation equations in Sec. IIIII and 
present our results in Sec. [IVJ 



II. BASIC IDEA OF CAUSAL HYDRODYNAMICS 



In this section, we review the idea of causal hydrodynamics using the charge diffusion example. Our 
discussion here is heuristic, but it serves a good starting point since the dispersion relation used in 
this paper in fact takes the same form as this simple example as shown in App. IA 31 Readers who are 
familiar with the idea of causal hydrodynamics may skip this section and may go to App. [S] directly 
for technical details. 

The basic set of equations is the conservation law and the constitutive equation (Fick's law for the 
charge diffusion): 

= , Ji = -Dd iP , (2.1) 

where D is the diffusion constant. These two equations lead to the diffusion equation: 

d oP - Ddfp = . (2.2) 

The diffusion equation is parabolic, which does not satisfy causality. In fact, the propagator of Eq. (|2.2p 
in (I + l)-dimensions is given by 

which has a small but nonvanishing value even outside the lightcone x > ct. 

In order to restore causality, one needs a hyperbolic equation such as the Klein-Gordon equation. 
The conservation equation must be true, so what is wrong is Fick's law. In fact, if dip = for t = 0, 
then Fick's law tells that the current vanishes immediately, i.e., Ji(t) = for t > 0. However, one 



Our computations do not show which modes are reliable to obtain t w . However, just before we submitted the first 
version of this paper, Ref. |31| appeared, which studies the similar problem as ours. They argue that the shear 
mode is unreliable to obtain t^, but the sound mode is reliable. Thus, we use our sound mode results for physical 
interpretations. 
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expects that the current should die away in reality. In order to incorporate this effect, Fick's law may 
be modified as 

TjdaJi + Ji = -Dd iP , (2.4) 

where tj is a new transport coefficient. In this case, one obtains Ji{i) = Ji(0)e _t / Tj . Thus, the 
parameter tj is the relaxation time for the charge current Ji. 

The modified law with the conservation equation leads to the telegrapher's equation: 

Tjd*p + d o p-Ddfp = , (2.5) 

which is a hyperbolic equation. The new term may become important at early time or for rapid 
evolution. Also, one can regard this as a higher order expansion of an effective theory. Hydrodynamics 
is just an effective theory with infinite number of parameters phenomenologically, so it is natural that 
new parameters arise. 

A propagating solution p cx e -™*+ l 9 z leads to the dispersion relation 

-tjw 2 -iw + Dq 2 = . (2.6) 

The wave-front velocity can be estimated by taking the q — > oo limit: Wf IO nt = y/D/rj. Thus, 
the equation is consistent with causality if Uf r0 nt < c (For large w. higher order terms may become 
important though). Let us also consider the opposite limit q — > 0. Then, the dispersion relation for 
the hydrodynamic pole (whose dispersion relation satisfies w(q) — > as q — * 0) is written by 

w = -iDq 2 ~iD 2 Tjq 4 + 0(q e ) . (2.7) 

(Only one solution is compatible with the low-energy limit.) 

Israel carried out a systematic analysis @ , but the resulting constitutive equations are still compli- 
cated. We restrict the case of linear perturbations and decouple each modes. Hydrodynamic modes 
are decomposed as follows: 

{longitudinal mode (diffusive mode) 
transverse mode 

{longitudinal mode (sound mode) 
transverse mode (shear mode) 
transverse traceless mode 

Not all modes have a hydrodynamic pole since such a pole arises due to a conservation law. The 
standard transport coefficients appear in the following modes: the charge diffusion constant D in the 
diffusive mode (as is clear from the above example), the shear viscosity r\ in the shear mode, the bulk 
viscosity £ and the speed of sound v s in the sound mode. 

In addition, Israel introduced 5 new transport coefficients: Three are relaxation times for the diffu- 
sive, shear, and sound mode, respectively {tj,t v ,th, respectively). The other two are the couplings 
among different modes (a coupling between the diffusive and the sound mode ao, and a coupling 
between the diffusive and the shear mode a±). At the end of the day, the dispersion relations for the 
diffusive and the shear mode just take the form \2. 7| ) for the telegrapher's equation. [See Eqs. (|A54[) 
and (|A59[) in App. IA 31 ] We determine these coefficients from the gravity computation below. 

III. GRAVITY COMPUTATIONS 

According to the standard AdS/CFT dictionary, the bulk gauge field acts as the source for the 
global R-chargc current on the dual field theory. (The It-charge is a global charge which presents in 
SYM; In this sense, it is an analog of the baryon number in QCD.) Similarly, the bulk gravitational 
perturbations act as the source for the stress-energy tensor on the dual theory. Thus, our aim is to 
solve these bulk field equations. 
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Let us consider the bulk perturbations of a p-brane which take the form 

A„(r) e- ,wt+iqz , h^ir) e- iwt+i * z , (3.1) 

where z := x p . The perturbations can be decomposed by the little group SO(p — 1) acting on 
x % (i = 1, • ■ • ,p — 1). The gauge field perturbations are decomposed as 

scalar mode (diffusive mode): Aq, A Zl A ri 
vector mode: Ai. 

Similarly, the gravitational perturbations are decomposed as 

scalar mode (sound mode): h u , h tz , h zz , h k , h rr , h tr , h zr , 
vector mode (shear mode): h t i, h Z i, h r i, 
tensor mode: hij — Sijh k k /(p — 1). 

Such a decomposition is essentially the same as hydrodynamics above. Each hydrodynamic mode 
couples to the corresponding bulk perturbation. 

Many authors solve such perturbation equations in various backgrounds. 2 Our aim is to get the 
subleading corrections by regarding causal hydrodynamics as an effective theory expansion in higher 
orders. In this paper, we consider the diffusive mode, the shear mode, and the sound mode. 



A. Backgrounds 

In this paper, we compute transport coefficients for the SAdS p +2 backgrounds. These backgrounds 
appear as the "near- horizon" limit of various branes. The p = 3 case corresponds to the D3-brane 
which is the M = 4 SYM. The p = 2 case corresponds to the Dl-brane or M2-brane in 11-dimensional 
supergravity. The p = 5 case corresponds to the D4-brane or M5-brane. The Dl and D4-branes are 
dual to the 2 and 5-dimensional SYM with 16 supercharges. 3 

The SAdSp+2 metric is given by 

dr 2 

ds 2 p+2 = f(-hdt 2 + dx%) + — , (3.2) 



where 



/ - , 0.3) 



ro 



The surface gravity is given by 



h=l-(f) ■ (3.4) 
P+ 1 r Q 



K = 



R 



2 



(3.5) 



Since some of the backgrounds can be interpreted as the Dp-brane, let us directly consider the 
Dp-brane for completeness. The Dp-metric consists of a (p + 2)-dimensional metric and warped S s ~ p . 
According to Ref. [13], the dimensional reduction of the metric into the (p + 2)-dimension gives 

„ 2(8- P ) f i „ „ , , i dr 2 

ds 2 p+2 = r^r- J Zp{-hdt 2 + dx 2 p ) + Z 1+ p— > , (3.6) 

g 2 s = r-TZ- 1 - 1 ? , (3.7) 



2 For example, at the lowest order in w and q, the perturbation equations have been first solved in Refs. |l8t fiotl for 
SAdSs and in Refs. [27J, [2l|j for SAdS4,7. See also Ref. pi for a recent application of SAdS4. 

3 For a recent discussion of this duality, see, e.g., Ref. |23H and references therein. 
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where 



The surface gravity is given by 



/ T 



-(7-p) 



z ={r) ' ^ 



ro 



7-p 



h = i-hr ■ ( 3 - 9 ) 



5-p 

2 i? 2 ^ 



(3.10) 



B. Field equations (diffusive and shear mode) 



Our computation closely follow Rcf. |24j. Let us start from the diffusive mode. This amounts to 
solve the Maxwell equation: 



v M (v=^n = o 



(3.11) 



where \/—g := \Z~-g~p+2/ '<?eff' ^he effective coupling g a R may be position-dependent. 

It is convenient to introduce a new radial coordinate u (u := rg/r for even p and u := r^/r 2 for odd 
p). We choose the gauge A u = and use a Fourier decomposition: 



Then, the Maxwell equation becomes 

g 00 wA' - gg zz A' z = , 

d u {V^g 00 g uu A' ) - V^g 00 g zz {wqA z + q 2 A ) = , 
du {V=gg zz g uu <) - V~gg m g zz (wqA + w 2 A z ) = o , 

where ' = d u . From Eqs. (|3 . 13[) and (|3.14[) . one gets a decoupled equation for A' : 

'd u (V^g 00 g uu A' o y 



d_ 

du 



-99" 



+ 



oo 



In the backgrounds we use, the equations of motion reduce to the following form: 



d_ 

du 



u^h—iu^A') 
du 



v "I it) 



A' Q = 



(3.12) 



(3.13) 
(3.14) 
(3.15) 



(3.16) 



(3.17) 



where h = 1 — u v ; a, (3, and v, are the constants which depend on the backgrounds; to and q are w 
and q normalized by surface gravity k [or temperature T := k/(2tc)]: 



to 



27tT' 



2ttT 



(3.18) 



Incorporating the "incoming wave" boundary condition at the horizon u = 1 and asymptotic form 
u = 0, 



A' =C(l-u 



F(i 



(3.19) 



The function F(u) is a regular function whose form can be obtained perturbatively as a double series 
in tt) and q 2 : 



F(u) =F Q + tvF 1 + q 2 Gi + ro 2 ^ + roq 2 #n + q 4 G 2 + ■ ■ 



(3.20) 



G 



Geometry 


shear mode 


diffusive mode 


Q 




V 


Q 


/J 


V 


SAdSp+2 (p: even) 


P 


-p 


P + 


1 


P-2 


2-P 


p+1 


SAdS p+ 2 (p: odd) 


P+1 

2 


P-1 
2 


P + 
2 


1 


P-1 
2 


3-P 
2 


p + 1 
2 


Dp (p: even) 


3 


p — 6 


7- 


P 


P-2 


-1 


7-P 


Dp (p: odd) 


2 


p — 5 
2 


7- 
2 


P 


P-1 
2 





7-P 
2 



TABLE I: The parameters a, (3, and i/ appeared in Eq. (|3.17[) . 



The dispersion relation is obtained by imposing Dirichlet boundary condition at it = 0. Such a disper- 
sion relation has been obtained for various theories at 0(tr,q 2 ). Our task is to compute corrections 
at O(ro 2 , roq 2 , q 4 ) to get the coefficients of causal hydrodynamics. 

For the shear mode, denote one of x l coordinates as x. We consider a metric perturbation of the 
form h tx ^ 0, h zx ^ with the other h^ v — 0. As explained in Ref. [24|, the equation for the shear 
mode reduces to the Maxwell equation. First, consider a fictitious Kaluza-Klcin compactification along 
the x-direction. Following the standard procedure of the Kaluza-Klein reduction, set Aq = (gzz)^ 1 ^^ 
and A z = (g zz )~ 1 h zx . Then, write the resulting action in terms of the Einstein metric. In the end, 
the perturbation equation becomes Eq. (|3.16[) with the replacement 

4-9 -> QzzV^ , (3.21) 

and the only differences arc the parameters a and (3 in Eq. (|3.17|) . 

Using the backgrounds in Sec. IIII A) one obtains a, (3, and 7 as shown in Table Q] It is easy to see 
that the perturbation equations for the Dl and the D4-brane are identical to those for SAdS4 and 
SAdSy, respectively. 



C. Field equations (sound mode) 

We closely follow Ref. [25j . The sound mode has 7 degrees of freedom. Out of these, there are 3 
gauge freedoms and 3 constraints, which leaves us a single degree of freedom. Our task is to constitute 
a master field which represents this degree of freedom and to obtain the master equation. 

First, let us compare our notations with those of Kodama and Ishibashi [25| . They take the SAdS p +2 
metric as 

dr 2 2M 
ds 2 = -f(r) dt2 + J^+ r2 d °l = 9ab(y) dy a dy b + r 2 (y) da 2 p , /(r) = -Xr 2 - — , (3.22) 

where (a, b) = (t, r). As a result, our notations and theirs are related as follows: 



Kodama and Ishibashi Ours 
A -1/R 2 

M r p Q +1 /2R 2 

f(r) f(r)h(r) 

da 2 dxp/R 2 
k Rq 
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where k and q represent wave numbers in their and our notations, respectively. 

They have written down the equations for 3 gauge-invariant variables X(r), Y(r), and Z{f) 
[Eqs. (2.24a)-(2.24d) in Ref. 



dX 
dr 
dY 
dr 



— =X 
dr 



r h dr 



1 d(fh) 



2(/ h) dr 



(X-Y) 



r 2 fh (fh) 2 

o 

in* 

z , 



z 



(fhy 



= a{r) X + /3{r) Y + 7 (r) Z , 
where we define a, f3, and 7 by 

1 /?'o\ 4 fr \P~ 3 



a(r) 
(3(r) 
7(r) 



2 2 , P 

w • r ■ H 



i? 



2 2 



q 2 r 2 h + 



P 



1 



1 



P - 1 



3-1 /r 
2 W 

r 



,p- 3 



(3.23) 

(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 
(3.29) 



Here, X, Y, and iT are Fourier-transformed as in Eq. (|3.12p . Equations (|3.23[) - p.25p are coupled first- 
order differential equations for 3 variables, but they reduce to a second-order differential equation for 
a single variable by a constraint (|3.26[) . which is the master equation. 

However, the master equation derived by Ref. [25| [Eq. (3.5) in their paper] is not particularly useful 
for our purpose. In order to solve the eigenvalue problem as a series in in and q, it is necessary that 
one can take the limit to, q — > not only for the perturbation equation but also for the boundary 
condition. Namely, suppose that the lowest-order solution as a series in to and q has some fall-off 
behavior as r — > 00. This fall-off behavior must coincide with the one for the full-order solution. 
Unfortunately, the master equation derived in Ref. [25[ does not satisfy this criterion, so one must use 
a new master field and obtain the master equation for such a field for which one can take the limit 
to, q — > 0. 

After some trial and error, we found that the following form of the master field $ is useful: 



#(«) 



h(s) 



1 



P 



1 



■2p 



p +1 X(s) + 



P 



1 



2p 



; p+1 Y(i 



where 

s := r /r . 

From Eqs. (|3.23[) - (|3.26[) . the equation for $ is schematically written as 

du 2 







B 1 (u)— +B {u)$ 
du 



For even p, 
Bx{u) 
B (u) 



2p{p + (2p + 3)u p+1 } to 2 + q 2 { (p - l)u 2 ^ +1 ) - {p 2 + 7p + 2)u p+1 - 2p 2 } 
uh [ 2pK> 2 + q 2 {(p- l)uP+ 1 - 2p} ] 
1 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(P+1) 2 

4 h 2 [2pw 2 + q 2 {(p-l)uP+ 1 -2p} ] 

x [8 P W 2 u 2p - 4(p + l)q 2 u 2p + 2p to 4 + q 4 {(p - l)u 2 ^ +1] - (3p - l)u p+1 + 2p) 

-to 2 q 2 {4p - (3p - l)u p+1 }] , (3.34) 
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where u := s. For odd p (= 2p' + 1), 

B ,x _ p{b + 2K , + 1 +^}^ 2 + q 2 {y^ , + 1 )-(^ 2 +6y + 3K'+ 1 -py} 



fl («) 



(p' + l) 2 1 



4 uh 2 {ptv 2 + q 2 (p' uP' +1 - p)} 
x [4p to V - 4(p' + l)q 2 u p + p tr 4 + q 4 ft (p - p' t/ +1 ) 

-» 2 q 2 {2p-(3p' + lK' +1 }] , (3.36) 

where u := s 2 . Finally, incorporating the boundary condition at the horizon u = 1 and asymptotic 
form u = 0, set 

$ = /i- 1 - lrt, / 2 .sf+ 1 F(w) . (3.37) 

For SAdSs, the master equation for the sound mode has been obtained by Kovtun and Starinets 
[2(| in the way consistent with the hydrodynamic limit [Eq. (4.35) in the paper]. The above master 
equation is identical to theirs. Also, for the Dp-brane, the master equation has been obtained by Mas 
and Tarrio [13] [Eq. (3.17) in the paper]. The above master equation for SAdS4 and SAdSy is identical 
to those for the Dl and the D4-brane, respectively. Our variable F is related to the variables by the 
others as follows: 

Kovtun and Starinets: Z 2 (u) = hT im/2 u 2 F{u) , (3.38) 
Mas and Tarrio: Z (u) = h-" v/2 F(u) , (3.39) 

where Z 2 and Zq are the variables used in the papers above. 

In general, the perturbation equations in the sound mode are harder to solve at the second order 
than the ones in the shear mode. Moreover, the generic SAdS p+ 2 case is harder than the SAdS.5 
case. In order to simplify our analysis, we employ the following method. First, anticipating the 
hydrodynamic dispersion relation, set 

W =doq + diq 2 +d 2 q 3 H , (3.40) 

and obtain the solution as a series in q: 

F(u) = F Q + qFi + q 2 ^ + • • • . (3.41) 

The constant di is obtained at each order by imposing Dirichlct boundary condition on the solution 
Fi. Then, we solve the equation for i^+i using fj (with the determined constant di). 

IV. RESULTS 

A. The shear mode and the diffusive mode 

The solutions F{u) are rather cumbersome expressions, so we do not write them explicitly. (The 
solution for the shear mode and the sound mode of the Af = 4 SYM are written in App. [S]) The 
Israel- Stewart theory has 5 new constants (tj, t„, rn, oto, oti)- For the backgounds with no R-charge, 
the gauge field and the metric perturbations decouple: This implies that ao = a\ = 0. In addition, 
for conformal theories rn = due to the vanishing £. [See Eq. (|A65[) .] Thus, the main interests are 
tj and 1V 

From the shear mode of the Af = 4 SYM, we get 

ft '=-4- zi ^ q4+o(q6) - 
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Comparing Eq. (|4.ip with the dispersion relation (|A54|) . we obtain the familiar result r)/s = l/(47r) 
and 

1 - In 2 . , 

The other results are summarized in Table [II] and Table IIIIl For the diffusive mode, one gets the 
diffusion constant D 

as well as rj. 

For the Af = 4 SYM, rj has never been obtained, but the result is obvious a priori from a result 
of Ref. [HI and our dispersion relation for the diffusive mode (|A58[) . Also, Ref. [Tt} has computed 
TV for the Af = 4 SYM from a somewhat different setting. They consider an expanding plasma and 
obtained t„ which is 3 times smaller than our result. We believe that this is due to a missing term 
in their constitutive equation. The term is negligible for a plasma near equilibrium, but it is not 
negligible for the expanding plasma. In fact, the discrepancy is gone once one adds the extra term in 
the constitutive equation (See Ref. [29[ for details.) 



B. The sound mode 



The sound mode is interesting in the sense that the relaxation time tv appears in this mode as well. 
[See Eq. (|A68|) .] Since rn = for conformal theories, one can deduce tv from this mode as well, but 
here one encounters a puzzle. From the sound mode of the Af = 4 SYM, we get 

q iq 2 3-21n2 3 4 ^ 

w = 7rT + ^r< +0(q) - (44) 

Comparing this with the dispersion relation (|A68|) . one obtains 

2 - In 2 . . 

T * = - w ' (45) 

which does not agree with the answer obtained from the shear mode (|4.2p . The other results are 
summarized in Table IIVI One can see similar discrepancies for the other SAdS p+ 2 backgrounds as 
well. 

We have not located the origin of the problem. But, first of all, the Israel-Stewart theory is not 
the unique formalism for causal hydrodynamics. One well-known alternative is the "divergence type 
theories" fiol [TTJ] . The discrepancies we found may imply that the gauge theory plasmas do not really 
fit into the framework of the Israel-Stewart theory. 4 

For the Dp-brane, the dispersion relation is identical to the one of the corresponding SAdS solution. 
However, the hydrodynamic interpretation is different partly due to the different spacetime dimen- 
sionality They are nonconformal, so tjj ^ and one cannot determine 7V and rn separately from the 
sound mode alone. For the Dl-branc, 

18 - (91n3 - V3ir) , . 

(No for the Dl-brane). For the D4-brane, using Eq. (|A67[) . we get 

„ + IS = 2 36 -( 91n3 + ^^0.29fm. (4.7) 
15 45ttT v ; 



4 One interpretation is proposed in Ref. [3lll . which appeared just before we submitted the first version of this paper. 
According to the paper, the shear mode is unreliable to obtain -rvr, and one should use the sound mode. Thus, we 
will use Table ITTT1 for physical interpretations hereafter. 
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Geometry 


Dispersion relation 




SAdS 4 (M2 & Dl) 

SAdS 5 (D3) 
SAdS 7 (M5 & D4) 


w= a," .3(91n3 + V3 7 r) q4 + 0(q6) 

tn = -iq 2 -i(ln2)q 4 + 0(q 6 ) 

Wq> 8(91n8-Var) , + 0(B) 
4 64 


9 In 3 + V3?r 
24ttT 

In 2 
2^rT 

9 In 3 - \/3tt 
24ttT 



TABLE II: The relaxation time r.j computed from the diffusive mode. 



Geometry 



Dispersion relation 



7v from shear mode 



SAdS 4 (M2) 
SAdS 5 (D3) 
SAdS 7 (M5 & D4) 



iq 2 . 9-91n3 + VStv 4 , 6 , 



iq 2 . 1 — In 2 4 6 . 
t» = -JL-i q 4 + 0(q 6 ) 



iq 2 , 18-9In3- v^tt 4 , n . 6s 
re = — o - * To q + 0(q ) 



9-(91n3-^7r) 
24ttT 

1 - In 2 



2ttT 

18- (91n3 + ^7r) 
24vrT 



TABLE III: The relaxation time tv computed from the shear mode. 



Geometry 



SAdS 4 (M2) 
SAdS 5 (D3) 
SAdS 7 (M5) 



Dispersion relation 



s/2 4 



iq 2 , 15 - 91n3 + V3tt_ 3 , ^,_ 4 , 



48^ 



-q J +0(q 4 



3-21n2 3 n . 4 , 

+ — q + 0(q ) 



q »q 

V3 3 6^ 



q 2iq 2 24 - 9 In 3 - \/3tt 3 _ . 4 , 

ro = 7I^ + t&s q + ° (q) 



Tvr from sound mode 



18- (91n3- V3tt) 



24ttT 
2 - In 2 



0.18 fm 



2ttT 

36- (91n3 + \/37r) 
24ttT 



0.21 fm 



0.27 fm 



TABLE IV: The relaxation time computed from the sound mode. Numerical values shown correspond to 
T" 1 = 1 fm. 



C. Discussion 



The theories we consider here are not QCD. Thus, it is important to ask if there is any universality 
or any generic behaviors just like rj/s. This is the reason why we consider various theories. From 
Table IIII1 there seems no obvious universality, but the numerical values of tv are similar among the 
theories we consider. 

To be more specific, get some numbers. First, recall that he ~ 197 McVfm and 197 MeV is not 
far from the QCD transition temperature T c . This means that the characteristic length scale at T c is 
T -1 ~ O(fm), and this is the typical value one would expect for relaxation times. In fact, the kinetic 
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theory predicts that 0] 

t " = T p = tS (48) 

for a 4-dimensional Boltzmann gas, where we used e+p = Ts and used the fact that the Af = 4 theory 
is conformal so that e = 3p. If one uses the AdS/CFT value of -q/s = l/(47r), = 3/(2ttT) ~ 0.5fm 
for T _1 = 1 fm. This value is not far from the our results in Tabic [TV] 

This seems consistent with what Israel and Stewart found. They obtained Eq. (|4.8[) by analyzing 
the Boltzmann equation. More precisely, they estimated P2 = T 7r/(2r/). 5 They found that P2 is not 
sensitive to the value of the cross section. This implies that @2 is more or less constant as we vary the 
coupling constant. Namely, r\ strongly depends on the coupling, and so does 7V, but T^/r] does not 
strongly depend on the coupling. (The entropy density s does not strongly depend on the coupling 
[HI, so it is irrelevant here.) And in fact, we found that the ratio T^jr] from the AdS/CFT duality is 
not far from the kinetic theory estimate. 
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Note added: While this paper is in preparation, a number of interesting papers appeared [30L l3ll. 
[Hj], which study the similar problem as ours (See also Ref. [33j|). In particular, Ref. [3l[ partly uses 
the same technique as ours. For SAdS5, our results of r T coincide with the results of Ref. [3l| both in 
the shear mode and in the sound mode. Also, Ref. [3l[ argues that the shear mode result is unreliable 
and one should use the sound mode to extract r^. Based on their observation, we have changed an 
early interpretation based only on the shear mode. See Ref. [13] for a review. 



APPENDIX A: ISRAEL-STEWART THEORY 

Sections I A 1 1 and I A 21 review the Israel-Stewart theory; Readers who are familiar with the formalism 
may go to Sec. I A 31 directly. 

1. Preliminaries 

Denote the number of spatial dimensions by d s . We consider the case of only one conserved charge 
p for simplicity, but the case of several charges is straightforward. The fundamental variables in 
hydrodynamics are the conserved current j 1 *, the energy- momentum tensor ', and the entropy 
current (which gives the direction of time) . 

a. Equilibrium 

In equilibrium, there is a special "fluid rest frame" defined by u£ q (u^ q = —1), in which there is no 
(spatial) flow. Thus, 

f = P< q , - e < q < q + Pcq (<r + < < q ) , s» = s cq < q . (Al) 



5 /?2 is a parameter they use to parametrize causal hydrodynamics [Eq. 1A2H 1 . However, /?2 appears in the dispersion 
relation only in the combination 2r?/32, which is tv. [Eq. | |A19| |1 
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The first law T eq ds cq = de — p cq dp tells that ,s oq is not an independent variable, but it is a function of 
e and p. Also, the temperature T cq and the chemical potential p cq arc defined by the first law as 

J_ = ds cq (e,p) Meg = ds oq (e,p) 

T cq de ' T cq ' dp ■ { ) 

The pressure p cq is not independent cither due to the Eulcr identity p cq = — e + T cq s cq + p cq p. It is 
convenient to rewrite the Eulcr identity in a covariant manner: 



*" = 7^-(Peq < ~ "cq,, - Mcq ^) . (A3) 



b. Near- equilibrium 

We closely follow Ref. @ but use slightly different conventions and notations. Consider a state of 
near-equilibrium whose deviation 6 from the equilibrium is small. In equilibrium, the entropy density 
s cq is a function of the charge p and the energy density e. We assume that the entropy current s M is 
a function of the currents j M and T M1/ even in the case of a noncquilibrium state: 

For a noncquilibrium state, various currents have spatial flows and they do not match in general. 
Thus, the notion of the "fluid rest frame" is ambiguous: a different current defines a different "fluid 
rest frame." There are two common choices for the "fluid rest frame" in the literature (The notations 
are defined below): 

1. The Eckart frame or Particle frame (N- frame): = in this frame. 

2. The Landau-Lifshitz frame or Energy frame (E-frame): fc^ = in this frame. 

Instead of choosing a particular frame, wc consider a general reference frame (u 2 = — 1) which is 
close to a fictitious rest frame of equilibrium thermodynamics, where 

| u» - < q | = 0(6) . (A4) 

Then, we derive the results so that they do not depend on a choice of ( "ffame-invariance" ) . Namely, 
we use neither the Eckart frame nor the Landau-Lifshitz frame. (However, we frequently comment 
the case of the Landau-Lifshitz frame since it is frequently used.) 

Given u M , one naturally defines a (d s + l)-decomposition of the spacetime tensor g^ v by the pro- 
jection operator h^ u (u): 

~ g »" + u" . (A5) 

Then, and T^ u are decomposed as 

j" = ptt"+# J (A6) 
= £/ 1 i"+2 k { l u v) + (p + n) W" v + . (A7) 

Here, the variables with "_L" represent the components which are orthogonal to u^: e.g., u^j^ = 0. 
(ir^ also satisfies u^ir^ = 0.) The quantities with "*" are the quantities defined by the functional 
form of the entropy density s eq (e, p) in the equilibrium. For example, use Eq. (|A2| for T(u) and fi(u). 
By construction, the "net flow of charge" jf_ (u) , the "energy flow" (u) , the trace and traceless part 
of the viscous stress H(u) and tt^ v (u) should satisfy 

|tf |, \k^\, |H|, \^\=0{5) . 
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We will make use of the "covariant time derivative" V u := u^V^ and the "covariant spatial deriva- 
tive" Z? M : e.g., D^v x = h ^ 'h u CT V ' p v x for a spatial vector v x . 
Also, one can write 

Q 

V^u,, = -Up a v + — h pu + cr M „ + uj^u , (A8) 

using 

:= V^u" = Wftfuvr , (A9) 
a" := V u u^ , (A10) 

:= (fy/ - V„u„ , (All) 

c«V := fy/ V p w CT . (A12) 

These quantities represent the expansion, acceleration, shear, and rotation of the reference frame, 
respectively. 

Now, the fundamental variables 7 71 , T' 1 ", and do not depend on a particular frame, but their 
(d s + l)-decompositions depend on a frame. Thus, let us check which variables are frame-invariant. 
Let us consider the following transformation of the reference frame: 

u" — it" := (l + Cl) 1/2 ^ + C^ , (A13) 
where u 2 = -1 , u„C±_ = , | CI | = C(<5) . (A14) 

One can check (See Appendix of Ref. Q) 

• The variations are 0(6 2 ): p, e, II, and n^ u . 6 

• The variations are 0(6): j x and k x . 

Thus, the former are the frame-invariant quantities up to O(S). The variables j x and fc^ are not, but 
they combine to give frame-invariant quantities up to 0(6): 

J» := f x --^W Xi (A15) 

IP* := u " + . (A16) 

e +p 

Note that J M consists partly of j x and partly of fc^. In the Landau-Lifshitz frame where k x = 0, 
_ an( j f/M — u m Using these variables, one can rewrite the currents as 

j» = P U" + J" , (A17) 

= (e + p) U y - ^ + pg» u + Tlh fiu + tH 4 " . (A18) 



2. Entropy and the second law of thermodynamics 

The constitutive equations are constructed so that the second law of thermodynamics V p s M > 
is guaranteed. Therefore, the form of the entropy current s M (as a function of and T M ") becomes 
important. 



Obviously, p, /i, and s are also frame-invariant as well since they are defined through p and e. 



14 



Assume that the Euler identity (|A3[) remains a good approximation even for noncquilibrium states. 
Thus, define Q^{u) by 

s"=: j (puV-UvT^-fij^-QV , (A19) 

where \Q»\=0(5). 

Fro small deviations, it suffices to retain only the 0(S 2 ) terms for Q^. We take the form 7 

=: + R" , (A20) 

where 

Q" = ^r(/? II 2 +/3iJ- J + ^TT^Tr^+^nj^ + ^TT^JA , (A21) 

R» := ^ [ 7r" A fc ± A + n/cV + — k± ■ k± \ ■ (A22) 

T(e+p)\ 2 J 

Note that i? M vanishes in the Landau-Lifshitz frame. Five constants a a and Pa appeared in Q^: (3 a 
represent the relaxation times and cxa represent the couplings among various modes as we will see 
below. 



a. "First order formalism" 

The divergence of Eq. (|A19[) is given by 

TrM" 4. TT / ft \ 

V = V M U u - J M W |J - V M Q M + 0(5 3 ) . (A23) 

Here, we have taken i? M into account (See Sec. I A 2 bj) . In order to ensure the second law of thermo- 
dynamics, the right-hand side of Eq. (|A23|) must be positive-definite: This is how the constitutive 
equations are derived. When one ignores Q M ("first order formalism"), V M s M > is guaranteed if the 
right-hand side is a sum of complete squares. Thus, introducing the transport coefficients, D, £, and 
77, we require 

J* = -Df (^y 1 D x (i) , (A24) 



9P, 

n = -ce, (A25) 

tt^ = -2?7£^, (A26) 

where 

/ h h pa \ 
9 := h^V„U v , SV := U(/Kf ~ — d J V p U a . (A27) 

The variables and E^„ are the frame-invariant expansion and shear, respectively. Then, 



= — r hf + — + - - V M Q" + 0(5') , (A28) 



Since we used frame-invariant quantities, it is clear that Q is frame-invariant up to 0(<5 2 ). It is convenient to construct 
Q in this way since V^s M = 0(<5 2 ). On the other hand, R M is not frame-invariant. The term _R M is chosen such that 
it cancels unphysical O(S) terms in V^s' 1 , which appear in the "first order formalism" below. 
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so V^s^ 1 > is ensured if D > 0, ( > 0, and r\ > 0, provided that (dfi/dp)f > 0. 
Note that the stress tensor is rewritten as 

T ij />//'■' s (-)//'•' -2 //>:'•' , 

using Eqs. (|A25|1 and (|A26j) . This is just the familiar form for the stress tensor. 



(A29) 



b. Second order formalism 

"The first order formalism" is not a closed form since one has to take into account the 0(5 2 ) terms 
i? M . 8 If one does not include this term, V M s M becomes frame-dependent; Moreover, hydrodynamic 
equations become an overdetermined system. On the other hand, if one includes this term as was 
done in Sec. IA 2 a[ one had better include all 0(S 2 ) terms in Q M . The second order formalism takes 
into account. In this case, the divergence of the entropy current gives 



Ju 



T 



T 



(A30) 



where Vjj := J7 M V M . Then, the constitutive equations for the second order formalism are given by 



J A = -d 



dp j T 

n = -c (e + ^Vt/n + aoV^j^ 

7r^ A = -2r? 



L> A ( 4 ) + /3i h x p VuJ p + a D x n + ati (D v tt vX + ir Xv a, 



h%h x a -^-p^ 



(3 2 VuK pa + ai D {P J^ 



(A31) 
(A32) 
(A33) 



3. Dispersion relations 

a. Assumptions and tensor decomposition 

In order to obtain the dispersion relations for causal hydrodynamics, we make a number of simpli- 
fying assumptions: 

• Linear perturbations: We consider linear perturbations from the thermal equilibrium. 

• Rest frame: We choose the fluid rest frame in equilibrium as the reference frame so that 

= , a„ u = w„„ = . (A34) 
Moreover, all thermodynamic quantities have no time-dependence and there are no flows, so 
J x = = , = u" , SU" u„ = . 

(The boldface letters represent background values and (W^ 1 := — U 11 .) 

• Rat (boundary) spacetime: We consider the flat (d s + l)-dimensional spacctime. Then, a M = 0. 



The Landau-Lifshitz frame is free of this problem since W = 0. 
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• "Decoupled ansatz": Let us take into account the bulk results in advance. In the text, we consider 
the backgrounds with no R-charge, so one can set p = jx = 0. Also, it holds (dp/dp) e = 
(d((iT~ 1 )/de)p = 0, and the gauge field and the metric perturbations decouple. This implies 
otQ = cti = 0. 

Introduce the time coordinate by u 11 = (dt) , and denote spatial indices as i,j, k, ■ ■ ■ . The linear 
perturbations are defined, e.g., by p =: p + 5p. Using Eqs. (|A17|) and (|Al8j) . one gets the conservation 
equations 



= d t Sp + pSQ + D i 6J i , 
= d t 5e+ (e + p) 5Q , 

= d t 5U l + — ^— ^ 

e + p 

Using Eqs. (|A31|) - (|A33|) . one gets the constitutive equations 

-l 



8J % 



-D 



dp 
dp 



TD l 5 



/3i d t 6J l + a D l SU + anDjSirV 



SU = -C [ 60 + /3b d t 5n + a Q D 3 5J J ] , 

h» h kl 



8ir lJ 



-2r? 



(A35) 
(A36) 

(A37) 



(A38) 

(A39) 
(A40) 



We henceforth use the "decoupled ansatz." 

In order to solve these equations (|A35|) - (|A40[) . decompose the spatial tensors as follows (A := DjD^): 



5 J 1 = D\J L + 4 , 
5U l = D l V L + U l T , 



St: 1 - 



A) 7r T + 2 J D ( Vi ) 



(A41) 
(A42) 

(A43) 



Here, J^, U^, Tt l T , and irji represent the transverse components which satisfy 



DiJip = 



DiU l T = , 



DiTT l T 



0, 



(A44) 
(A45) 



Note that tcl = since Stt 1 ^ is traceless. 



b. Vector modes 

The vector modes consist of U^, J^, and tt^. The equations for these fields are given by 

= d t U* T + , (A46) 

e + p 

= tj d t J l T + 4 , (A47) 

= T7r d t n % T + tt 1 t + n U l T , (A48) 

where the relaxation times are defined by 

T^:=2r]02, (A49) 

t~j := D p\ . (A50) 



dp 



T 
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Eqs. (|A46[) and (|A48|) combine to give 

= t„ d 2 t -K l T + d t TT l T - ^-r . (A51) 
The Fourier-Laplace transform of the equation 

f(w, q) := /°° d£ / dz e -™*+^ /(t, z ) , (A52) 



Jo J 

gives 

,2 



= -t^w 2 ~iw+ T ^- . (A53) 



Here, we used the Euler identity e + p = T s for zero chemical potential. 
Expanding Eq. (|A53[) for small w and q 2 , one obtains 



w = -i D n q 2 - i D 2 t v q 4 + 0(q 6 ) , (A54) 



where D v := rj/(Ts) 



c. Scalar mode (diffusive mode) 



There are two scalar modes: the diffusive mode and the sound mode. We discuss them separately 
since they decouple due to the "decoupled ansatz." 
The diffusive mode consists of Sp and Jl- One gets 

= d t Sp + AJ L , (A55) 
TjdtJL + JL+DSp, (A56) 



which combine to give 

= TjdfSp + d t 5p- D A5p . (A57) 

This takes the same form as the shear mode equation (|A51|) . so one can immediately write the 
dispersion relation: 

= -tj w 2 - i w + D q 2 , (A58) 

or 

w = -iDq 2 -iD 2 T.,q 4 + 0(q e ) . (A59) 

d. Scalar mode (sound mode) 

The sound mode consists of Se, Ul, SH, and ttt- One gets 

= d t 5e+ (e + p) MJ L , (A60) 

= 8 t U L + V s 5e + m + ( 1 : 1 ^ A * T , (A61) 
e + p 

= rn d t Sn + 5R + C AU L , (A62) 

= 7v d t TT T + ttt + 2 n U L , (A63) 
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where 

2 d P 
m ■= (fa ■ 

The Fourier-Laplace transformation of Eqs. (|A60[) - (|A63[) gives 
= (1 — irn w)(l — i 7V w)(w 2 — v 2 q 2 ) 



iwq 
e + p 



Thus, the dispersion relation for the hydrodynamic pole is 



i ( d s — 1 C\ 

w = ±v s q — — 7] + -)q 

e+p \ d s 2/ 



±- 



1 



d, - 1 



1 - 1/4 



V +C « s T"n - 



1 - 1/da 



'2i> s (e +p) [ (i s \ " e+p 
In particular, for conformal theories where £ = Tn = 0, 

, 1 - 1/d. 2 1 1 - 1/4 ( 2 

to = ±w s q - i ; r? g ± — — ») 2t) a r T 



e+p 



2 v s e + p 



e+p 4(e + p) 

i - l/d 



e+p 



7? g J + 0(g 4 ) . 



(A64) 
(A65) 



(A66) 



(A67) 



(A68) 



e. Tensor mode 



The tensor mode consists only of irj! , and its equation is given by 

TV 9t7Ty + 7Ty = . 

Thus, the dispersion relation for the tensor mode is 

i 



w 



, 



but this is unreliable since it is inconsistent with the hydrodynamic limit. 



APPENDIX B: PERTURB ATIVE SOLUTIONS FOR AT = 4 SYM 



(A69) 
(A70) 



Here we give explicit expressions for the perturbative solutions for the Af = 4 SYM. Integration 
constants are fixed by requiring the solutions to be regular at the horizon. For the shear mode, 9 

C 



F = C, F l = iC[ --+ln(l + u 

u 



2u 



(Bl) 



F 2 = C 
H n = - 
G 2 = lc 



-l + ln2 In 2 



iC 



+ ^ In(« - 1) - J In(u + 1) hi{(l-u) 2 (u + 1)} - ^ Li 2 



u + 1 



1 - lln(u+ 1) ) +ln \-\ 
u \ 2 v ' I \u+l 



1 / u 
- + In 

u \ u + 1 



(B2) 
(B3) 
(B4) 



Note added in v6: In previous versions, we erroneously listed the solution for the diffusive mode. However, F\ as the 
diffusive mode solution is still incorrect and should be F\ = Cln(l + u)/(2u). 
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where Li2(u) is a polylogarithm. 
For the sound mode, 

F = \ Cl (l + izM) , (B5) 
and the Dirichlct boundary condition at u = determines do = =bl/v3- Using these, we get 

and the boundary condition gives di = — i/3. Finally, using these lower order results, we get 

F 2 = 3_C[-ln 2 (u + l)u 2 + 21n21n(w- l)u 2 -21n(l - u) ln(u + l)u 2 + 4 In 2 ln(u + l)u 2 

u + 1 



-2Li 2 



8u± 24^3 d 2 -81n(u + l) +81n2- 12 



(B7) 



and the boundary condition gives d 2 =±(3-21n2)/W3. 
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